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ABSTRACT 

Cyclotron decay and absorption rates have been well studied in the literature, focusing primarily on 
spectral, angular and polarization dependence. Astrophysical applications usually do not require reten- 
tion of information on the electron spin state, and these are normally averaged in obtaining the requisite 
rates. In magnetic helds, higher order quantum processes such as Compton scattering become resonant 
at the cyclotron frequency and its harmonics, with the resonances being formally divergent. Such di- 
vergences are usually eliminated by accounting for the finite lifetimes of excited Landau states. This 
practice requires the use of spin-dependent cyclotron rates in order to obtain accurate determinations 
of process rates very near cyclotronic resonances, the phase space domain most relevant for certain ap- 
plications to pulsar models. This paper develops previous results in the literature to obtain compact 
analytic expressions for cyclotron decay rates/widths in terms of a series of Legendre functions of the 
second kind; these expressions can be expediently used in astrophysical models. The rates are derived 
using two popular eigenstate formalisms, namely that due to Sokolov and Ternov, and that due to John- 
son and Lippmann. These constitute two sets of eigenfunctions of the Dirac equation that diagonalizc 
different operators, and accordingly yield different spin-dependent cyclotron rates. This paper illustrates 
the attractive Lorentz transformation characteristics of the Sokolov and Ternov formulation, which is 
another reason why it is preferable when electron spin information must be explicitly retained. 

Subject headings: radiation mechanisms: non-thermal — magnetic fields — relativity — stars: neutron 
— pulsars: general — gamma rays: theory 



1. INTRODUCTION 

The number of astrophysical sources thought to possess sur- 
face magnetic fields above the quantum critical field strength 
of B CI = m 2 e c i /{eh) w 4.413 x 10 13 G has been steadily grow- 
ing in recent years. The five presently known Soft Gamma- Ray 
Repeaters (SGRs) and six or seven Anomalous X-Ray Pulsars 
(AXPs) are believed to be magnetars, neutron stars having sur- 
face fields in the range 10 14 - 10 15 G (Duncan & Thompson 
1992). In addition, the Parkes Multibeam survey (Manchester 
et al. 2001) has discovered at least six new radio pulsars that 
have surface fields near or above critical, and several with fields 
comparable to those of some of the magnetars. The radiation 
processes that determine the observed X-ray and gamma-ray 
spectra in these sources are operating in the extreme relativistic 
and quantum regimes and thus require treatment that is ac- 
curate in such environments. One process that is particularly 
important in the source emission models is resonant Compton 
scattering, in which electrons scatter photons at the cyclotron 
resonance with a cross section much larger than at continuum 
energies. Relativistic electrons can blue-shift low energy pho- 
tons into the resonance, upscattering the photons into a high- 
energy continuum (Daugherty & Harding 1989, Dermer 1990). 



The quantum clectrodynamical (QED) cross section for 
Compton scattering at the cyclotron fundamental (Herold 1979) 
and higher harmonics (Daugherty & Harding 1986, Bussard, 
Alexander & Meszaros 1986), accurate in arbitrarily high mag- 
netic fields, has been known for some time. However, these 
derivations do not treat the natural line widths that render the 
cross section finite at the cyclotron resonances. In order to use 
such rates in spectral calculations in astrophysical models, the 
line widths that originate from the lifetimes of the excited Lan- 
dau states must be included in the cross section (Wasserman & 
Salpeter 1980). The width of the nth resonance is equal to the 
cyclotron decay rate from that state (Pavlov et al. 1991), and 
the prescription for incorporating the widths in the QED cross 
section in the high-field regime has been discussed by Harding & 
Daugherty (1991) and Graziani (1993). Including resonant line 
widths in the scattering cross section necessarily requires spin- 
dependent decay rates, which appear in an infinite sum over 
Landau state n and spin of the intermediate virtual states, even 
in the case of ground state-to-ground state scattering in the fun- 
damental. Different spin states have different decay rates, and 
thus different resonant energy denominators. Since the elec- 
trons and positrons in intermediate states have non-zero mo- 
mentum parallel to the magnetic field, it is important to use 
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basis states that yield a spin dependence that is Lorentz invari- 
ant, i.e. boosts along the magnetic field do not lead to a mixing 
between the spin states. 

As is true for all quantum processes, the spin-dependent rates 
and cross sections depend on the choice of electron wavefunc- 
tions in a uniform magnetic field. Historically, several choices of 
wavefunctions have been used in calculations of the scattering 
cross section and cyclotron decay rates. The two most widely 
used wavefunctions are those of Johnson & Lippman (1949) and 
Sokolov & Ternov (1968). The Johnson & Lippman (JL) wave- 
functions are derived in Cartesian coordinates and are eigen- 
states of the kinetic momentum operator. The Sokolov & Ternov 
(ST) wavefunctions, specifically their "transverse polarization" 
states, are derived in cylindrical coordinates and are eigenfunc- 
tons of the magnetic moment operator. Given the different spin 
dependence of the ST and JL eigenstates, one must use caution 
in making the appropriate choice when treating spin-dependent 
processes. Herold, Ruder & Wunner (1982) and Melrose and 
Parle (1983) have noted that the ST eigenstates have desir- 
able properties that the JL states do not possess, such as being 
eigenfunctions of the Hamiltonian including radiation correc- 
tions, having symmetry between positron and electron states, 
and diagonalization of the self-energy shift operator. As found 
by Graziani (1993), the ST wavefunctions also diagonalize the 
Landau-Dirac operator, and are the physically correct choices 
for spin-dependent treatments and for incorporating widths in 
the scattering cross section. Although the spin-averaged ST 
and JL cyclotron decay rates are equal, the spin-dependent de- 
cay rates are not, except in the special case in which the initial 
momentum of the electron parallel to the field vanishes. 

Perhaps the most fundamental argument against use of the 
JL wavefunctions is that radiative corrections cause excited JL 
Landau levels to become unstable to spin-flip transition within 
Landau states on timescales comparable to the timescale for de- 
cay to a lower Landau n state. These radiative corrections to 
the magnetic moment break the spin degeneracy of both the ST 
and JL excited Landau states (Herold, Ruder & Wunner 1982), 
causing spin-dependent energy level shifts, and implying a de- 
pendence of the rates for spin-flip transitions within a Landau 
state on the choice of wavefunctions. The rates of such spin-flip 
transitions within split ST Landau states have been evaluated 
(Geprags et al. 1994; see also Parle 1987), and are found to 
be of order of af (£?/_B C r) 4 for B <C B CI , which is negligible 
compared to decay rates between Landau states. The existence 
of such relatively long-lived states is a premise of the S-matrix 
formalism and is also essential to any astrophysical calculations 
involving spin-dependent transitions between Landau states. In 
contrast, it is anticipated that since the JL eigenstates do not 
diagonalize the self-interaction Hamiltonian including radiation 
corrections (i.e. the mass operator), the associated mixing in- 
curred in S-matrix evaluations will render spin-flip transition 
rates with fixed n comparable to ordinary cyclotronic rates with 
changes in n , a situation that is unphysical. 

In this paper we discuss the Lorentz tranformation character- 
istics of cyclotron decay rates for both ST and JL formulations, 
and derive simplified expressions for the decay rates from an 
arbritrary excited Landau state to the ground state. We show 
that the ST eigenstates preserve separability of the spin de- 
pendence under Lorentz boosts along the local magnetic field, 
a desirable property that does not extend to the JL formal- 
ism, for which such Lorentz boosts mix the corresponding spin 
states. In the ST formulation, by taking advantage of its spin- 
state preserving characteristics, our analytic simplifications can 



be elegantly applied to the spin-dependent decay rates; in the 
JL formalism the simplifications can be compactly applied only 
to the spin-averaged rates. The resulting expressions replace 
integrals over emergent photon angle with series of Legendre 
functions of the second kind, which correspond to sums of ele- 
mentary functions, easily yielding simple asymptotic forms. 

The expressions derived here should have wide applicability 
to modeling cyclotron emission, Compton scattering and other 
QED processes in super-critical fields. For such magnetar-type 
fields, resonant scattering takes place primarily at the funda- 
mental, since the cyclotron energy exceeds 1 MeV. Although 
the resonant scattering line widths formally involve infinite sums 
over Landau states, in the case of the fundamental resonance the 
sum is dominated by the n — 1 state, whose width is equal to the 
n — > cyclotron decay rate. For cyclotron scattering in higher 
harmonics I > 0, the intermediate sums in the vertex functions 
have the largest contributions from the n = I + 1 state. How- 
ever, for B 3> -B C r j cyclotron transitions to the ground Landau 
state dominate (Sokolov, Zhukovskii & Nikitina 1973, White 
1974, Harding & Preece 1987), so that cyclotron decay rates for 
n —>■ transitions treated in this paper should be good approxi- 
mations to the widths of excited states, a circumstance pertinent 
to astrophysical models of magnetars. 

2. CYCLOTRON RATES: SOKOLOV AND TERNOV FORMALISM 

This presentation of spin-dependent cyclotron rates appropri- 
ately focuses first on the formulation generated using the trans- 
verse polarization eigenstates of the Dirac equation derived by 
Sokolov & Ternov (1968) for electrons in a uniform magnetic 
field B. Explicit forms for them can also be found in Herold, 
Ruder & Wunner (1982) and Harding & Preece (1987), whose 
exposition on cyclotron emission essentially forms the basis of 
results developed here. Herold et al. observed that these states 
diagonalize the operator that describes the electronic self-energy 
shift in an external magnetic field. Another attractive feature 
of the Sokolov & Ternov eigenstates is that they possess charge 
conjugation symmetries, i.e. between electron (positive energy) 
and positron (negative energy): see Eq. (12) of Herold, Ruder 
& Wunner (1982). An additional asset of these wavefunctions 
that is enunciated here is that they yield cyclotron rates whose 
spin dependence is effectively separable from Lorentz transfor- 
mations along the field: i.e. such boosts do not mix spin states. 
This is a very useful characteristic that is not present in the 
Johnson & Lippmann formalism addressed in Section [3] below. 

Herold, Ruder & Wunner (1982) obtained general expressions 
in their Eq. (17) for the spin-dependent cyclotron rate for tran- 
sitions between arbitrary Landau levels, but for the case of zero 
initial momentum p parallel to the field. Latal (1986) indepen- 
dently obtained similar results for cyclotron transitions to the 
ground state, but retained arbitrary initial momenta p along B; 
his Eq. (22) forms the starting point for the exposition here. De- 
noting £ = ±1 as the spin quantum number of the initial elec- 
tron, ground state transitions are characterized by a single final 
spin £' = — 1 . The energy level quantum numbers are initially 
n , and n' = after the transition. Latal's spin-dependent total 
rates can be written in the form 
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where B is the magnetic field strength, expressed dimension- 
lessly hereafter in units of the quantum critical field B CI = 
mlc 3 /(eh) « 4.413 x 10 13 Gauss, and the spin-averaged, total 
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cyclotron transition rate for n — > n' 
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Here A = h/m e c is the Compton wavelength of the electron over 
2n , and other quantities in this equation are defined as follows. 
For initial and final dimensionless momenta (i.e., in units of 
m e c), p and q respectively, parallel to the field, the electron's 
initial energy E n and final energy Eq , both dimensionless (i.e., 
in units of m e c 2 ), are 



E n = y/l + 2nB + p 2 and E = y/T 



(3) 



This convention of using dimensionless energies and momenta 
will be adopted throughout the paper. 

The integration variable /i is the angle cosine of the emit- 
ted photon with respect the magnetic field direction, i.e. \x = 
cos9 = k.B/\k\ \B\ for photon wavenumber vector k . No inte- 
gration by parts has been performed when obtaining Eq. J5J, so 
that the cyclotron rate, differential in photon angles, dT n o/dfi , 
again averaged over initial electron spins, corresponds directly 
to the integrand of Eq. (J2J : 



r 



dp 



(4) 



Corresponding spin-dependent differential cyclotron rates can 
be deduced from Eq. Q in similar fashion. 

The remaining variable, £ , emerges naturally from energy- 
momentum conservation in the interaction. Only momen- 
tum parallel to the field is conserved, since the system is not 
translationally-invariant orthogonal to B. Hence cyclotron emis- 
sion satisfies two such conservation relations: 



En 



P = 



E +uj 
q + uJ/d 



(5) 



for \x = cos 9 . Here u> is the dimensionless photon energy, hav- 
ing been scaled by m e c 2 . The elimination of q from the system 
in Eq. (J5J leads to the identities 
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1 + £ E n - pfi ' 
c^(l-/i 2 ) = (l-£)(E n -pp) 
that are simultaneously satisfied, with 
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2nB(l 
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(6) 



(7) 



being the n' — > specialization of the definition in Eq. (lib) of 
Latal (1986). The angle integration in Eq. J5J for the cyclotron 
rate is non-trivial, due largely to the complicated dependence of 
£ on fi , and the presence of the exponential. 

This completes the definitions relevant to Eq. However, 
we also note that the identities in Eq. (JSJ) combine to yield 



^ 2 (1-M 2 ) 
2B 



1-g 



(8) 



which is precisely the argument of the exponential in Eq. 
This exp(— k) factor results from the n' = specialization 
of associated Laguerre functions that are formed from the spa- 
tial integration (Fourier transform) of the eigenfunctions. Since 
co sin 9 is Lorentz invariant under boosts along B, so also is n . 
Furthermore, forming E n — pp using Eq. JSJ), and then elimi- 
nating lj using the second identity in Eq. |JB}, quickly leads to 
the equivalence 

E - qp 



£ = 



E n - PV 



(9) 



When weighted by u) top and bottom, this is a ratio of invariant 
products of four-momenta associated with the interaction; i.e. 
£ is also a Lorentz invariant for boosts along the field, which 
can be inferred directly from Eq. JHJl. 

2.1. Lorentz Transformation Properties 

There is no need to compute Eq. (0 as it stands for arbitrary 
p . Instead, one can appeal to a simple Lorentz transformation 
protocol. Consider boosts along the field between the inertial 
frame where the initial electron possesses momentum p along 
the field, and that frame where the initial parallel momentum 
is zero. In this latter "rest" frame let fio be the photon an- 
gle cosine with respect to B, po = be the initial parallel e~ 
momentum, and the initial electron energy be denoted 



(10) 



e n = VT+ 2nB . 

Then, for a dimensionless boost velocity (3 = v/c and Lorentz 
factor 7 = l/yl — f3 2 along B, the initial electron quantities 
in Eq. 05J satisfy 



E n — 
P = 



l(po ~ f3s n ) 



7^n ; 
-7 far. 



(11) 



from which (3 = —p/E n = —pj y/l + 2nB + p 2 is established. 
The photon angle cosine satisfies the aberration formula 



/i = cos f = 



1 - @Po 



(12) 



This latter relation identifies a suitable change of variables 
for the angle integration, for which d/i = dp,o/[j(l — /3/io)] 2 ■ 
Other transformation identities include u) sin 9 = lu sin 9q and 
(1 — n 2 )/(E n — p[i) 2 = (1 — /io)/ e n I from which the invariance 
in form of £ is established via the definition in Eq. J7J). It then 
follows, after a modicum of algebra, that this pure change of 
variables leads to an alternative form for the rate in Eq. ©: 



=; CKfC 
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dpo (i - 



(i+o 



n— 1 / i\ / 

— (4 - cxp( -n j 



i-e 



(13) 



(14) 



where, now £ is obtained from Eq. Q using the correspon- 
dences /i — > po , p — > and E n — > e n ■ 

Accordingly, this identifies attractive Lorentz transformation 
behavior, with the rate reduced (i.e. lifetime dilated) by just 
the Lorentz factor 7 of the boost. Moreover, such boosts along 
the field keep the spin states "separated," i.e. there is no im- 
plied mixing of states incurred by such transformations. This 
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inherent simplicity is an appealing characteristic of the Sokolov 
& Ternov states, and was identified by Graziani (1993); it is 
not exhibited in the Johnson & Lippman formalism explored 
in Sec. 01 Clearly, this extraction of p > cases via a sim- 
ple modification factor outside the integral is expedient for the 
subsequent analytic developments. 

2.2. Analytic Developments 

While the integral expressions for the cyclotron rates can be 
routinely evaluated numerically, they can also be represented by 
compact analytic series in terms of elementary functions that 
are readily amenable to computation. The integration variable 
fio for I n (B) in Eq (|14|) is not the most convenient; a more 
expedient choice for the purposes of analytic reduction is 



K 

n 



1-g 




(15) 



The integration is even in fj,Q , and since 1 / e n < £ < 1 , the 
integration range maps over to 



< 



< 



1 



1 



(16) 



After a modest amount of algebra, the change of variables then 
leads to the form 



In{B) 



,-n<t> in-1 



\J {4>n ~ 4>) (1/07) 



1 



(17) 

This integral can be evaluated in terms of Appell functions, 
the degenerate, two-dimensional hypergeometric functions, us- 
ing Eq. 3.385 of Gradshteyn & Ryzhik (1980). Such a step does 
not facilitate evaluation, since Appell function expositions usu- 
ally develop double, infinite power series representations (e.g. 
see Burchnall & Chaundy 1941; Exton 1976). 

Analytic progress is quickly made via Taylor series expansion 
of the exponential e - ™^ around = 0. The order of inte- 
gration and the infinite summation can then be interchanged 
because the integration is uniformly convergent on the interval 
< <j) < <p n , since 4> n is strictly less than unity. Another change 
of variables <f> = e~ f , with the definition 



z n — coshirj 
then establishes 

In{B) 

with 



1 



-) s 1 ' 7^7 ■ (1S1 



E 

fc=0 



V2 Jt„ \J cosh t — cosh t r , 



1 — e ' cosh t n 



(19) 



(20) 



Using the identity 8.715.2 of Gradshteyn & Ryzhik (1980) gives 
an integral representation of the Legendre function Q y {z n ) of 
the second kind: 



Qu{z n ) 



V2 J tn \J cosh t — cosh t r< 



Zn > 1 



(21) 



These special functions are finite sums of elementary logarith- 
mic and polynomial functions of z n (e.g. see Abramowitz & 
Stegun 1965). Note that the Q v {z) are generalizable to associ- 
ated Legendre functions Q£(z) , defined in 8.703 of Gradshteyn 
& Ryzhik (1980), from which Q v (z) = Q°(z) . It follows from 
Eqs. iJUl and © that 



J n-\-k (^n) — Qn-\-k— 1 \%n ) Z n Qn-\-k (^n ) 



(22) 



and relevant properties of the Q v and J„ functions are listed in 
the Appendix. Computational issues for the J m are discussed 
there also, and an expedient approach for integer m = n + k is 
to compute the J m using the recurrence relation in Eq. l|A4fl . 
The final result of these analytic reductions is 
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k=Q 



(-n) k 
k\ 



Jn-\-k (^n) 



(23) 



which expresses the spin-averaged cyclotron rate as a compar- 
atively simple series of elementary functions. For low values 
of n , this series is rapidly convergent for all z n > 1 , typi- 
cally requiring 3-4 terms to achieve 0.01% accuracy. Such a 
rate of convergence renders the series evaluation computation- 
ally much more efficient than a numerical integration of Eq. (J2J 
by Simpson's rule or a quadrature technique. The appearance of 
the Legendre functions of the second kind, Q v [z) , in the series 
expansion for the rate is a result of the cylindrical symmetry 
imposed by the presence of the external magnetic field. 

Spin-Averaged Cyclotron Rates 
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Fig. 1. — The spin -averaged cyclotron decay rates r„o (solid curves) 
resulting from Eq. 1231 . or equivalently Eq. J2J, for different values of n , as 
labelled. The rates are scaled by the typical lifetime to = afc/A, and are 
plotted as functions of the magnetic field strength B , expressed in units of 
Bcr = 4.413 X 10 13 Gauss. The dotted lines represe nt th e power-law asymp- 
totic dependences (i) in the li mit nB < 1 in Eq. 1261 . for n = 3 , and (ii) 
in the nB S> 1 limit in Eq. 1281 . specifically for n = 30 . The dashed 
curve is the logarithm of the ratio of sp in-d ependent rates for n = 1 in the 
Sokolov & Ternov formalism, i.e., Eq. 1241 . Results are presented for zero 
electron momentum initially along B, i.e., p = or 7 = 1. 
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Evaluations of Eq. ll'.'il) are presented in Fig. for different 
values of n and for 7 = 1, demonstrating the well-known rapid 
increase of the rate with B when B <; 1 . The dependence 
on B is relatively weak for supercritical holds. The Figure also 
clearly exhibits the property dT n o/dn < 0, which can, with 
extensive effort, be deduced from the integral form for the spin- 
averaged rate in Eq. {T3J|. This decline of r n0 with increasing 
n is evident in the asymptotic approximations for nB <C 1 and 
nB 1 derived shortly. Fig. also exhibits the logarithm of 
the ratio 

r-J = VI + 2nB + 1 

r^" VI + 2nB - 1 V ' 



of the spin-dependent rates for n — 1 ; higher n cases display 
similar behavior. This curve illustrates how the rate is approx- 
imately independent of initial spin when nB 1 , whereas 
£ = — 1 initial states lead to much faster transition rates (by of 
order (nB)^ 1 ) when nB <C 1 , reflecting the well-known dom- 
inance of non-spin-fhp cyclotron transitions by non-relativistic 
electrons (e.g. see Melrose & Zheleznyakov 1981). 

2.2.1. Asymptotic Limits 

Now consider the two pertinent asymptotic limits for the spin- 
averaged rate in Eq. When nB < 1 , (\> n « nB/2 < 1 
and z n > 1 . Using the hrst identity in Eq. (JXTJ, Eq. 8.771.2 of 
Gradshteyn & Ryzhik (1980) can then be invoked to obtain the 
leading order dependence of J„(z„) , noting the specialization 
2-Fi(a, P, 7, 0) = 1 of the common hypergeometric function. 
The result, for general arguments z , is 



(n + l)r(n)r(l/2) 
2r(n + 3/2) 



(2*)" 



z > 1 



(25) 



Here, T(x) represents the Gamma function. This is the lead- 
ing order (fc = 0) term in the series evaluation of I n (B) , and 
hence this result is readily checked by taking the limit 

of Eq. I|17|) ■ The doubling formula for the Gamma function, in 
Eq. 8.335.1 of Gradshteyn & Ryzhik (1980), can now be used to 
yield an asymptotic form 



r„ « ^(2n a ) n l2nV()\ B 



(n + 1)! 



n+1 



nB < 1 



(26) 



which is just Eq. (31d) of Latal (1986). This form illustrates 
the strong dependence of the cyclotron rates on both n and B . 
A sample n = 3 case of the limit in Eq. (|26|) is displayed in 
Fig. ^ This limit can only be compared with the well-studied 
classical cyclotron limit in a restricted fashion, since it is still 
an essentially quantum result. Moreover, it incorporates the in- 
trinsically relativistic effect of spin-orbit coupling. Observe that 
Eq. I|26|) . in conjunction with Eq. Q), reproduces both the non- 
spin-fhp (£ = — 1 ) and spin-flip ( £ = 1 ) ST rates encapsulated 
in Eqs. (18) and (21), respectively, of Melrose & Russell (2002), 
or equivalently, the n' = specialization of Eq. (18) of Herold, 
Ruder & Wunner (1982). 

Classical formulations of cyclotron emissivities (e.g. see 
Eq. (6.19) of Bekcfi 1966) yield rates that depend on compo- 
nents (3±c of the electron velocity perpendicular to the field: 
r„ oc B0 1 ™ • Here the classical rate 

r„ effectively represents the nB <C 1 limit of a sum over 
various n — > n' transitions, with n' < n . By invoking the 
correspondence 2nB — > p\ = j3j_ , it is quickly established 
that r n oc B n+1 , the same dependence as in Eq. (|26|l . Yet 
in general, the classical r„ does not equal the quantum r n0 , 



since here n' > contributions are omitted from considera- 
tion. However, in the particular case of n = 1 , there is only 
a single n' < n final state, namely n' = , and with the pre- 
scription f3j_ — > 2B , Eq. (6.19) of Bekefi (1966) specializes to 
yield T% = 4ctfC B 2 / (3X) , which is precisely the value obtained 
in Eq. I|26|l . Derivation of the nB 3> 1 , ultra-quantum limit is 
slightly more involved. In this situation, <f> n — ► 1 and z n — > 1 
so that the entire k series in Eq. I|23|) is retained. The key useful 
result is derived in the Appendix, in Eq. (|A5|1 , which establishes 
J n+ k{z) — ► l/(n + k) as z — ► 1 + . Including the next order log- 
arithmic contribution then leads to an approximate form for 
I n (B) in Eq. (|A6|) . so that it then follows that 




l e~ n log. 2nB 



rw 



2nB 



nB > 1 



(27) 

Here, 7(11, x) is the incomplete Gamma function in 8.350.1 of 
Gradshteyn & Ryzhik (1980). Note that the same result can 
be derived directly from the <j> n — ► 1 limit of Eq. fTTjl using 
the integral definition of the incomplete Gamma function. The 
weak \f~B dependence of the rate on B in this limit is now 
apparent, and is evinced in Fig. ^ by the crowding of curves for 
supercritical fields. In the limit of n 1 , further specializa- 
tion is possible using the limit 7(71, n)/T(n) — > 1/2 as n — > 00 , 
derived in the Appendix, to yield the approximation 



ate / B 

InO « — V8^ 



n > 1 , B > 1 



(28) 



This ultra-quantum result corresponds to the limit given in 
Eq. (32a) of Latal (1986), though the numerical coefficient of 
Latal's does not agree exactly with the 1/V8 factor in Eq. (|28|l . 
the difference being around 10%. However, the expressions for 
n — > quantum cyclotron transition rates derived in Eq. (9) of 
Sokolov, Zhukovskii & Nikitina (1973) and Eq. (2.56) of White 
(1974) reduce exactly to Eq. (|28|l in the limit B ^> 1 , confirming 
the asymptotic analysis here. A sample n — 30 case of this limit 
is displayed in Fig.^ the precision of which is around 2% relative 
to the exact result, and better than Eq. (32a) of Latal (1986). 
The first order correction to this asymptotic limit that would 
follow from the logarithmic term in Eq. I|27() is 0{\og e n/ yfn) 
relative to the leading order contribution in Eq. I|28l) . 

2.2.2. Differential Photon Spectra 

In the p — frame, the spin-averaged differential cyclotron 
rates dT n o/duo can be immediately extracted from Eq. Ill3() 
using the prescription implied by Eq. JIJ. However, given 
that Eqs. JJjJ and Q define a one-to-one correspondence be- 
tween the angle cosine n Q and the energy u of the emitted 
cyclotron photon, the differential rates can suitably be pre- 
sented as differential photon spectra, dT n $/duj. The photon 
energy is ui — (e n — l/e„)/(l + £) , so that the kinematic bound 
1/^n < £ < 1 translates to 



nB 



Vl + 2nB 



ui- < oj < ui + = Vl + 2nB - 1 



(29) 



The change of variables is then defined by two branches, 

Mo = ±77 V( £ n - 1 - w)(e n + 1 - u) , (30) 
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with the factors inside the square root reflecting the kine- 
matics associated with both the cyclotron radiation and the 
one-photon pair creation interactions. This has a Jacobian 
du /cLu — — £(1 + £,) £ n/[{ £ n — l)MoH , leading to the spin- 
averaged differential spectrum 



is the form for the cyclotron power for general p > , with 



duj 



a.fC 1 



A ~B^T(n + l) Vc^ 
(e„ — uj)nB — lu 

7^ 



= exp —n — 
l-w)(e„ + l-w) V w - 



(31) 



This spectral rate is summed over the polarizations of the ra- 
diated photons. Note that a factor of 2 has been included to 
account for the two branches of the relation between u and 
lu . Observe also that this spectrum can be obtained directly 
from Eq. ((13(1 and the form for I n {B) in Eq. (|17fl by using the 
correspondence (f> = lo/lj— — 1 . 

Key characteristics of the differential spectrum include its 
broadening and shifting to higher energies as nB increases, 
behavior that is largely governed by the kinematic result in 
Eq. ((291) . An exponential decline is present in n = 1 cases until 
a sharp kinematic peaking at lu ~ lu + = e n — 1 overtakes it. 
In n> 1 cases, this exponential contribution is dominated by 
the power-law factor (lu/lu^ — l) n , leading to a steeply- rising 
spectrum. When nB ^> 1 , nearly all of the radiative power in 
the spectrum is confined to energies cu ~ \/2nB . In the oppo- 
site domain, nB <C 1 , the spectra are extremely narrow, and 
little insight is gained by attempting an involved comparison 
with the classical limit; most classical results are rendered for 
sums over a multitude of transitions n — *■ n' . 



K n {B) 



o / duo (1 - g) 



i + 1 ) n+2 r u 



u exp i- n T+! 



( 34 ) 

where £ is given by Eq. 1(15(1 . This expression for the power 
is identical to that in Eq. (28) of Latal (1986), who illustrated 
the magnetic field and n dependence in his Fig. 1. The spin- 
dependent powers P^q are, of course, just given by additional 
multiplicative factors: P^ Q = (1 — (/e n )P n o for £ = ±1 . 

Observe that the 7 dependence hasjiow disappeared from the 
analysis, highlighting the fact that P n o is a Lorentz invariant 
for boosts along B. This follows from the exact compensation 
of time dilation and photon rcdshifting effects when integrating 
over all angles; such an attractive invariance property for the to- 
tal power clearly does not extend to its differential counterpart 
dPno/dfi . 

The integral for K n (B) can be manipulated in similar fash- 
ion to that for I n (B) . Since the only difference is an extra 
2/(1 + £) factor, which translates to a (1 + cf>) factor in the 
light of Eq. (USJ), it follows that K n (B) can be worked into 
the form of Eq. 1)17(1 . but with the extra (1 + 0) factor. The 
subsequent developments in Sec. l2.2l rapidlv lead to the identity 



K n (B) = ™' I Jn+k{Zn) + Jn + fc + l(Zn)} 



fe=0 



(35) 



2.2.3. Cyclotron Power 

To conclude the formalism based on Sokolov & Ternov states, 
expressions for the cyclotron energy loss rate, or power, _P„o — 
dE„o/dt can be expeditiously presented. The development of 
such rates parallels that for r„o , with the inclusion of an extra 
(dimcnsionless) weighting factor to to the differential decay rate 
dr n Q/du . Hence, one can quickly write down 



P 1l0 = m e c / uj ^77— da 



l d/i 



amieC n n B 



lu d/i 

(1-0 



(32) 



ra-l 



as the spin-averaged power, using Eq. J5J. This can be conve- 
niently transformed to the integration variable fio , representing 
the p — frame, using an analysis similar to that presented 
in Sec. 12.11 The only change is the manipulation of the ad- 
ditional lu = 2nB/ (1 + £)/(E n — p/i) factor inside the integra- 
tion, for which the Lorentz transformation relations in Eqs. (|ll|l 
and l(T2|l quickly yield l/(E n — pa) — 7(1 — [3a )/e n . Since 
dT n ()/daQ is even in uq , it follows that the term proportional 
to uq dT n o/daQ in the transformed integrand of Eq. I|32|l does 
not contribute. Therefore, comparison with Eqs. (|13fl and (|14H 
then indicates that 



P 



nO 



a{m e c n 



n+l 



B 2 



71 



IT! 



K n {B) 



(33) 



for insertion into Eq. ((33(1 . Given this form, asymptotic lim- 
its can be routinely obtained by replicating the approaches 
of Sec. |2~2~T1 When « 1 , it is evident that K n {B) w 
I n {B) , since then J n+ i(z n ) <§; J n (z n ) , as can be inferred for 
z n ^> 1 fro m E q. 1(25(1 . Furthermore, for nB ^S> 1 , the re- 
sult in Eq. ((A5() quickly yields a leading order contribution 
K n (B) « j(n, n)/n n + -f(n + 1, n)/n n+1 . Then the identity 
j(n + 1, n) = nj(n, n) — n n e~ n and Eq. ((A8() can be used to 
derive the limiting form as n — > 00 . The results are 



PnO 



2 n n 2n+l (n + 1) ! Bn+ 2 nB j 

aim e c J (2n + ly. 

■j , n, B > 1 . 



(36) 

While the nB <C 1 case here replicates that inferred from 
Eq. (31) of Latal (1986), the nB > 1 result differs from the 
numerical result in Eq. (32a) of Latal (1986) by around 10%. 

The mean cyclotron photon energy (lu) is just the ratio of 
P n o to r„o , so that one quickly deduces from a comparison 
of Eq. (HP with Eqs. and @ that (lu) w nB when 

nB < 1, as expected from classical cyclotron theory, and that 
(lu) w \/2nB when nB 1 . Note that Latal (1986) misses the 
V2 factor (see his Eq. (32b) for the nB ^> 1 domain) due to the 
imprecise nature of his estimates of the pertinent numerical fac- 
tors in P n o and r„o . These behaviors of (lu) are reflected by 
the differential spectral form in Eq. H31|) , where for nB <§; 1 the 
spectrum is essentially a delta function pinned at the cyclotron 
harmonic energy nB , and where for n ^S> 1 and B > 1 the 
spectrum is somewhat broad but markedly asymmetric, skewed 
strongly towards lu w lu + w V2nB . 
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3. THE JOHNSON AND LIPPMANN FORMULATION 

A complementary formulation of the cyclotron problem was 
provided by Daugherty & Ventura (1978, hereafter DV78) us- 
ing wavefunctions derived by Johnson & Lippmann (1949). This 
is essentially the foremost presentation in the literature of cy- 
clotron radiation calculations that specifically uses the John- 
son & Lippmann eigenstates, hereafter referred to as JL states. 
These states are eigenvalues simultaneously of the Dirac Hamil- 
tonian and the kinetic momentum operator 7? = p+eA/c , where 

A is the vector potential, and therefore are not eigenstates of 
any physically meaningful spin operator (see Melrose & Parle 
1983). The relative lack of symmetry in the JL states between 
the electron and positron wavefunctions leads to more compli- 
cated Lorentz transformation characteristics for the cyclotron 
rates for boosts along the magnetic field, as shall become evi- 
dent below. The spin-averaged rate should (and does), of course, 
reduce to the Sokolov & Ternov formalism result in Eq. J2J. 

3.1. Reduction to a Latal-Type Form 

For direct comparison with the Sokolov & Ternov formalism 
results above, it is necessary to generalize the cyclotron emission 
results in the Appendix of DV78 to treat arbitrary n — > tran- 
sitions, specifically their Eqs. (A2) and (A3). This is a routine 
S-matrix calculation in quantum electrodynamics, and essen- 
tially parallels the cyclotron absorption calculations presented 
in Sec. II of DV78. The reader is therefore referred to DV78 for 
relevant background and some definitional material, and here 
only a detailed outline of the derivation is provided. Wherever 
possible, the same notation as used in the ST formulation above 
is adopted. The starting point is Eq. (Al) of DV78, which can 
be integrated over photon solid angles dil = 2ir d(cos 9) to yield 
the spin-dependent cyclotron rate 



1 nO 



2tt 



QffC 



/ d(cos 6) dcu $<r 

=-L, II J 



(37) 



where dn/dcu = 4itL 3 uj 2 /(2irX) 3 is the density of photon states 
in a box of size L ; this size establishes the normalization of the 
wavefunctions. Note that dn/du) is dimensionless because of 
the adopted convention of dimensionless photon energies ui . 

A sum over the a =_L, |j photon polarization states has been 
effected, with the standard convention for the labelling of the 
photon polarizations being adopted: |j refers to the state with 
the photon's electric field vector parallel to the plane containing 
the magnetic field and the photon's momentum vector, while _L 
denotes the photon's electric field vector being normal to this 
plane. Also, 



1 J_ 

U[B cT 



f L dq 




db 


Sfi 


1 A 2nJ 




2tt 



(38) 



is the relevant integration over the final parallel momentum q 
of the electron and the coordinate b defining the mean location 
of the final electron wavefunction perpendicular to B. Here T 
is the normalizing time for the QED perturbation calculation. 
The central contribution to the integrand comes from 



Sfi 



I a { X 3 2tt 
2^ y/E n Eo 



8{E n -Eo-to) V fi 



(39) 



the S-matrix element for cyclotron emission that is analogous 



to Eq. (14) of DV78 that describes cyclotron absorption. Here 
Vfi = [ d 3 xe rk - 3 ul{q, b, x)M a u^(p, a, x) (40) 



is the Fourier transform of the spatial portion of the interaction 
amplitude, basically the vertex function of Melrose and Parle 
(1983), but with the temporal phase factors already extracted. 

The initial electron wavefunction is ti„ (p, a, x) , and can pos- 
sess either Q = ±1 spin designation (note that this spin notation 
differs from the s = ±1 notation adopted by DV78); the final 

electron wavefunction, vl(q, b, x) , is an n = , C — — 1 ground 
state, so that an explicit spin label is suppressed. Expressions 
for these spatial parts of the Johnson & Lippmann eigenfunc- 
tions of the Dirac equation are given Eq. (10) of DV78; these 
do not exhibit the same charge conjugation symmetry that the 
Sokolov & Ternov "transverse polarization" states do. 

In Eq. 140fl . M.„ is one of the two "polarization matrices" in 
Eq. (15) of DV78, representing contracted products of polariza- 
tion vectors and the Dirac 7 matrices. Since the Un and itj 
functions possess the dimensions of (length) ~ 3 / 2 , it follows that 
Vfi and consequently Sfi and <& a are dimensionless. Observe 
also that the temporal integration over the wavefunction prod- 
ucts has already been performed leading to the 6 function in 
Eq. (gnj) that expresses energy conservation. Finally, note that 
a defines the mean location of the initial electron wavefunction 
perpendicular to B, i.e. is a counterpart to b. 

The squaring of the matrix elements in Eq. (|38(l is performed 
using the approach of, for example, Bjorken & Drell (1964) for 
handling the S function. The result is 



Of 2ttA 3 

2uj ~LT 



S(E n -Eq-jS) 
E n E 



Vfi 



(41) 



in parallel with the developments leading to Eq. (19) of DV78. 
Integrations analogous to those in Eq. (|40|l have been performed 
analytically in DV78 for the cyclotron absorption case, result- 
ing in the appearance of associated Laguerre functions that are 
characteristic of magnetized QED calculations involving either 
free electrons or electron propagators. Here, since specialization 
to the n' = case is made, these functions reduce to simple 
combinations of power-laws and exponentials, viz., 



^2 w 



K 

n! 



2 • 2 
ui sin 1 

2B~ 



(42) 



in the notation of Eq. (18) of DV78. As these spatial integra- 
tions embody translational invariance along the field, they yield 
a S function expressing conservation of momentum parallel to 
B, i.e. establishing p = q + uj cos 9 , and completing the conser- 
vation laws in Eq. JSJ. 

A simple crossing symmetry manipulation of Eq. (16) of DV78 
facilitates the evaluation of Eq. 140(1. and also Eq. I|38|l . which 
can be compared closely with Eq. (19) of DV78. Using the 
identity in Eq. I|42|) . and reverting to the fi = cos 9 notation, 
a modest amount of algebra then yields relatively compact ex- 
pressions for the spin-dependent cyclotron emission rate in the 
Johnson & Lippmann formalism: 



< _ 



(E n — Eq) 



CtfC / , 

4A E n (E - qf i) (n-l)\ 

(£„ + i)(£ + i){A;j + Ai} . 



(43) 
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The different photon polarization and electron spin combina- 
tions yield the following contributions: 



A 



A7 1 



q 



p 



E Q + 1 E n + 1 



En + 1 



2nB 



g p \ oj{1 - p?) 

E + 1 E n + 1 J 2nB 



(44) 

2 



for the || polarization, and 



a: 1 = 



'7 



P 



1 E n + 1 



(45) 



for photons of _L polarization. It is now almost trivial to demon- 
strate that for the n = 1 case, Eq. (|43(l in concert with Eqs. (I44|l 
and (|45|l reduces to Eqs. (A2) and (A3) of Daugherty & Ventura 
(1978), when £ = +1 and £ = — 1 , respectively. Accordingly, 
Eqs. I|43H45[) serve as a generalization of the cyclotron emission 
rates presented in DV78 to arbitrary n — + transitions. 

The isolation of the photon polarization contributions in this 
development is an attribute that is absent from the exposition 
in the Appendix of Daugherty & Ventura (1978). Such an isola- 
tion expedites the algebraic developments that reduce this form 
into expressions that much more closely resemble those in the 
Sokolov & Ternov formulation, i.e. Eqs. and ©. The ana- 
lytic reduction is routine, but lengthy, involving the expedient 
use of several identities. Among these are Eqs. (007 which can 
be used to quickly derive 



E n ii - p = E p - q , £,{E n - pfj,) = E - qp 



(46) 



the latter of which is just a rearrangement of Eq. ©. For the 
( = +1 case, it is just a short path to yield the equivalences 



A 



En 



E n + 1 



En 



1 



u {E + l)(J, - q 

E n + 1 En + 1 



(47) 



and for £ = — 1 , a moderate amount of work to yield 



V 



2nB 



nB 

E n +1JE + 1 E n + 1 



p 



(48) 



With this assembly of tools, and various algebraic manipula- 
tions, one is lead to the following final forms for the spin- 
dependent cyclotron rates in the Johnson & Lippmann formu- 
lation: 

C ate n" B f 1 e~ K dfi (1 - Q"" 1 
L »° ~ A (n-l)\ E n }_ x (E n ~prf(l+d) n+1 C ' 

(49) 

where Eq. JHJ) gives k in terms of £ , and the spin-dependent 
factors in the integrand are 



X+i 



2nB 



2nB 



+ 1) \E n -pn 
2(l + 2nB- |) -X+i 



- E n (1 - 



(50) 



These are clearly forms resembling the integral that Latal (1986) 
derived, though indicating a slightly more complicated differen- 
tial distribution dr n0 /dfj, than for the ST case. The average of 
the factors is obviously 



X 



X+i + A— 1 



1 + 2nB 



(51) 



from which it is evident that the spin-averaged cyclotron rate 
for this Johnson & Lippmann formulation is identical to that 
of the Sokolov & Ternov formalism, i.e. Eqs. Q and ©, as 
should be the situation. Such degeneracy between formalisms 
also prevails for the power P n n radiated when averaged over 
the initial electron spins. 

Note that while the spin-dependent rates are generally differ- 
ent for the two sets of wavefunctions, in the limit p — > , it is 
clear from Eq. (|50|l that X^ —►(]. — £/s n )X so that the two for- 
malisms become degenerate in this limit. This is a consequence 
of the mathematical identity of the Sokolov & Ternov and John- 
son & Lippmann eigenstates in the special case of p = . 

3.2. Lorentz Boost Characteristics and Analytics 

As with the analysis in Sec. [21 it is instructive to discern the 
Lorentz transformation behavior of the Johnson & Lippmann 
rates. Such a development closely parallels that in Sec. 12.11 
again writing the relevant integrals in terms of variables in the 
frame in which p — . With the identities E n — ^e n and 
l/(E n — pfj,) = 7(1 — Ppo)/s n , it can be quickly established that 



X 



+1 



2nB 



1 e n I 



X 



2nB^f) 



(52) 



which can be used also for X_i 
cyclotron power considerations of Sec 



2 X - X +i . As for the 
2.2.31 since the spin- 



averaged integrand or dT n o / dfio is even in /zo , the term propor- 
tional to po in Eq. (|52|l contributes exactly zero to the integral 
in Eq. 149JI . This simplification clearly also applies to X-i . It 
then follows that an alternative form for the rate in Eq. (|49|) is 



OfC 1 



B 



1)! ISn 



£n(j£n + 1) 



In{B) 



(53) 



This form conveniently avails itself of the analytic developments 
of Sec. 12.21 so that no further reduction of integrals is necessary: 
the evaluation of I n (B) as a series of Legendre functions of the 
second kind in Eq. I|19|l is immediately applicable. Such a sim- 
plification cannot be applied to the spin-dependent cyclotron 
powers, a point that is addressed below. 

The Lorentz-transformed expression for the cyclotron rates in 
Eq. (|53|l is attractively simple, but is not as elegant as its equiv- 
alent for the Sokolov &: Ternov formulation, as embodied in 
Eqs. and ljl3|l . Here, boosts along the field do not introduce 
simple time dilation factors of 1/7 , essentially amounting to a 
mixing of JL spin states induced by such transformations. This 
less than ideal characteristic is a consequence of the fact that 
the JL wavefunctions are eigenstates of a somewhat convoluted 
spin operator that is not symmetric under charge conjugation 
(e.g., see Melrose & Parle 1983), nor is it manifestly covariant 
under boosts along the field. In contrast, the ST "transverse po- 
larization" wavefunctions are eigenstates of the component \i z 
along B of the magnetic moment operator pi = a — ijx {p+e A} 



9 



(here a and 7 are matrix vectors for the Dirac algebra), and 
such eigenstates exhibit both appropriate symmetry between 
positrons and electrons, and simple transformation properties 
for boosts parallel to the field. 

The sensitivity of the spin-dependent cyclotron rates, and in 
particular their Lorentz transformation characteristics, to the 
choice of spin eigenfunctions is naturally expected, since the 
inherently relativistic effect of spin-orbit coupling is explicitly 
incorporated in the Dirac equation in an external field, thereby 
emphasizing the interpretative importance of the choice of wave- 
functions. The greater elegance of the results in Sec. |3 is yet 
another argument in favor of usage of Sokolov & Ternov states 
for the cyclotron problem, and therefore for other QED pro- 
cesses where spin-dependent cyclotron rates are required. 

Note that the differential spectrum for spin-dependent con- 
siderations can be routinely obtained in a manner similar to the 
derivation of Eq. H31[) . observing that X can be replaced by 
X + i in Eq. (|52p. or X_i = 2X — X + i , as desired, and then 
using the substitution £ = (e n — l/s n )/u> — 1. 



JL to ST Cyclotron Ratio 




0.5 



-1 1 

Log 10 [nB] 



Fig. 2. — The ratios of the spin-dependent cyclotron decay rates result- 
ing fro m Jo hnson & Lippmann and Sokolov & Ternov formalisms, as given 
m Eq. JPJ. The rates are plotted as functions of nB , with B expressed 
in units of B CI = 4.413 X 10 13 Gauss, for different fixed values of 7 , as 
labelled, corresponding to parallel momenta p = —■yf3e n . The upper three 
curves correspond to f = +1 , while the lower three are for f = — 1 . The 
dotted line represents 7 = 1, where the ratio is unity for both f = ±1 . 

The ratios of the JL rates in Eq. (|53|1 to their ST counterparts 
are plotted in Fig. [21 being specifically given by 



1 nO 


JL 


1 nO 


ST 



C(7- 



/ 









(54) 



In general, this ratio is less than unity for £ = — 1 and greater 
than unity for £ — 1 ; it is also clearly a monotonically increas- 
ing function of 7 for £ = +1 , and decreases monotonically with 
7 for £ = — 1 . This ratio approaches unity for either 7 — > 1 
(dotted line) or for 2nB ^> 1 , and in the case of £ = — 1 also 
for 2nB « 1. In the limit 2nB <C 1 , the ratio approaches 



the maximum departure from unity in 
1 domain, as 7 — > 00 , for £ = — 1 the 
(4 + 2\/2)/(4 + 3V2) w 0.83 and occurs 



27/(7+1) when C = 1 ■ 
this case. In the 2nB ~ 
minimum in the ratio is 
for nB = 1 + V2 . 

The significant deviations from unity underline the impor- 
tance of the choice of wavefunctions when retaining spin de- 
pendence. This becomes a profound issue for resonant Comp- 
ton scattering problems (e.g. see Graziani 1993), where, as 
discussed above, the spin-dependent cyclotron lifetime is re- 
quired to render the intermediate electron states metastable and 
thereby truncate the divergent resonance at the cyclotron fre- 
quency. It is also pertinent to handling resonances in magnetic 
Coulomb scattering (e.g. Langer 1981), two-photon pair cre- 
ation in strong fields (e.g. Kozlenkov & Mitrofanov 1986), mag- 
netic photon splitting above pair creation threshold (discussed 
briefly in Baring 2000), and other general external field prob- 
lems involving electrons either in free states or as propagators 
(Graziani, Harding & Sina, 1995). 

To conclude this section, a brief exposition of the spin- 
dependent powers is offered. The spin-averaged power is, of 
course, that in Eq. (|33|l . To retain spin dependence in the JL 
formalism, one proceeds by weighting the integrand of Eq. I|49|) 
with a u> — 2nB/ (1 + £)/ (E n — pfi) factor. Guided by the anal- 
ysis in Sec. \2.2^ terms that are odd in fiQ contribute zero, and 
the by now routine manipulations yield the form 



P ( 



a.(m e c 



A 



1 - 



B 

C(7 + £n) 



(55) 



1) 



A'„( B ) + <|gf| i „(B, 



where K n (B) is the integral in Eq. (|34() . or equivalently, the se- 
ries in Eq. I|35l) . and the L n {B) term emerges from a /1q term in 
the integrand resulting from the product loXq , and corresponds 
to 

(-n) ' Jn+k(z n ) 



Ln(B) = ]T 



fe=0 



fc! 



fe + 1 



(56) 



Eq. I|55|) is clearly not as elegant as the corresponding ST re- 
sult, namely P^ = (1 — (/s n )P n o, but does reduce to this 
simpler form in the limit (3 — > , as expected. Furthermore, 
averaging Eq. over spins £ = ±1 then yields exactly P n0 , 
so that degeneracy between the JL and ST formalisms is real- 
ized when spin information is not explicitly retained. Note that 
asymptotic limits of P„ for small and large nB can be readily 
obtained but are not particularly elucidating for the purposes 
of this exposition. 

4. CONCLUSION 

This paper presents useful analytic developments of spin- 
dependent cyclotron decay rates/ widths for general n — ► tran- 
sitions between Landau states, expressing these in terms of vari- 
ous series of Legendre functions Q n {z) of the second kind. The 
rates and the radiated cyclotron powers are derived using two 
popular eigenstate formalisms appropriate for the Dirac equa- 
tion in a uniform magnetic field, namely that due to Sokolov & 
Ternov (1968), and that due to Johnson & Lippmann (1949). 
The resulting expressions are Eq. (|23() for the cyclotron rate and 
Eq. I|33|) in conjunction with Eq. (|35|l for the Sokolov & Ternov 
formulation, and Eq. ijSlfl) for the rate and Eq. l(53|) for the power 
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for the Johnson & Lippmann case. These compact forms can 
be expediently used in astrophysical models, since they gener- 
ally amount to rapidly convergent series of elementary functions. 
These spin-dependent results are particularly important for ad- 
dressing resonant cyclotronic divergences that appear in various 
higher-order QED mechanisms in external magnetic fields. The 
Johnson & Lippmann developments extend the work of Daugh- 
erty & Ventura (1978) to arbitrary n — ► transitions, while 
the series development permits a more accurate determination 
of nB 3> 1 asymptotics of the Sokolov & Ternov rates than 
that offered in Latal (1986). 

In addition, the Lorentz transformation characteristics of the 
rates are derived for boosts along the magnetic field, yielding 
a simple time dilation multiplicative factor in the ST case, but 
a slightly more complicated and less ideal dependence on boost 
Lorentz factor for the JL analysis, where spin-state mixing is 
implied by such boosts. Such comparitive complexity of the 
JL formalism is more pronounced for the radiated power. This 



comparison again underlines the superior attributes of Sokolov 
& Ternov eigenstates for use in evaluating rates and cross sec- 
tions of magnetized processes in quantum electrodynamics. 
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APPENDIX 

PROPERTIES INVOLVING THE ASSOCIATED LEGEND RE FUNCTIONS OF THE SECOND KIND, J V (Z) 

This Appendix provides several identities germane to the developments in the text, principally in relation to the Legendre functions 
of the second kind, Q v (z) , which form the basis of the series expansions for the cyclotron rates. They can be defined in terms of 
hypergeometric functions (e.g. see Eq. 8.820.2 of Gradshteyn & Ryzhik 1980), or alternatively the integral representation in Eq. (|21|l . 
For integer indices v = n , the specialization of interest here, they are simple combinations of logarithmic and polynomial functions 
(e.g. see Abramowitz & Stegun 1965), that lead to similar simplicity for the J v (z) using Eq. (|22l . Observe that Q v {z) = Q°„{z) , 
where the associated Legendre functions Q„(z) arc defined in 8.703 of Gradshteyn & Ryzhik (1980). Note also the alternative 
identities 



Ju{z) 



)„{t)dt 



z > 1 



(Al) 



embodied in the results 8.734.3 and 8.752.5 in Gradshteyn & Ryzhik (1980). Explicit functional forms for some of the Q n (z) are 
given in Abramowitz & Stegun (1965) and Eq. 8.827 of Gradshteyn & Ryzhik (1980), from which the following identities for the 
first few J n (z) can be formed: 



,h(z) 
Hz) 



= -i 



i 



log 



1 



3z 2 



z- 1 
2 



Zz[z z - 1) 



log 



z+l 



(A2) 



13z 
6 



5z J 



(5z 2 -l)(z 2 -l) 



log 



z + l 



An effective way to derive expressions for the J n {z) in terms of elementary functions, at least for low n , is to use the Rodrigues 
formula in 8.836.1 of Gradshteyn & Ryzhik (1980), and then invoke Eq. (|22p. An alternative approach is to use identities 8.831.2 
and 8.831.3 of Gradshteyn & Ryzhik (1980), to arrive at a relation between Q n (z) and a sum of Legendre polynomials P s (z) and 
a logarithmic term. This quickly yields 



Jn(z) 



Qn-l(z) - zQ n (z) 



(A3) 



P n . X {z)-zP n {z)\ log ^±l\-J2^Pk-l(z)P n -l-k(z) + zJ2^Pk-l^)Pn-k(z) , 



fe=l 



k=l 



a form that clearly displays the combined logarithmic and polynomial character of these functions, and which is useful for asymptotic 
results, explored just below. 

There are various possible approaches to numerical calculations of the J n functions for integer indices. Options include direct use 
of Eq. i|A3|l , combined with an efficient algorithm for computing the Legendre polynomials for arguments z > 1 , such as invoking 
recurrence relations (e.g. see Eq. 8.5.3 of Abramowitz & Stegun 1965), which are numerically accurate for upward recurrence in 
the polynomial index n. Alternatively, direct recurrence of the Q n (z) functions, which satisfy the same relations as the P n (z) , 
is viable, though this is only stable for downward iterations in n when z > 1 . This approach therefore requires evaluation of 
the P n (z) for two particular initial values n, perhaps via the hypergeometric function series using the identity in Eq. 8.1.5 of 
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Abramowitz & Stegun (1965). However, probably the most expedient numerical technique is to use the recurrence relation for the 
J n {z) functions directly: 

n(ri + 1) J n+1 (z) - n(2n + 1) z J n (z) + (n 2 - 1) J„_i(z) = . (A4) 

This identity can be deduced after a modicum of algebra from Eq. 8.5.3 of Abramowitz & Stegun (1965). Note that strong 
cancellation often arises between the last two terms of this equation, so that for large n and z > 1, the values of J n (z) decline 
exponentially with n according to Eq. (|25J) . Starting with initial values prescribed by Eq. (|A2|) . this technique is stable to upward 
recurrence in n for any z > 1 , the range of arguments appropriate for the cyclotron problem. This quickly yields suitably precise 
values (numerically tested for n < 100 ) for a range of indices required in the series summations that appear in the cyclotron rates 
and powers, for example in Eqs. ill'MI) and (|35|l . The use of this recurrence becomes superfluous typically for n > 30 , since such 
indices are generally required only in the classical cyclotron limit when z n > 1, and the asymptotic limit in Eq. I(25|) is then the 
preferred computational tool. 

Now focusing on asymptotics issues, two limiting values of J n (z) are required for the analysis in Sec. 12.21 The z » 1 limit is 
routinely derived, and is listed in Eq. (|25|) . Since z > 1 , the other limiting domain is z — > 1 . This is a somewhat more involved 
case. By inspection of Eq. i|A3l) . as z — > 1 , the logarithmic terms approach zero since P n (z) — > 1 + n(n + l)(z — l)/2 + 0(z — l) 2 . 
The Legendre polynomials thus simplify to approximately unity, to leading order, and it follows that in the neighborhood of z = 1 , 



Jn{z) = ~ - {z-l) log 



+ 0(z-l) 



(A5) 



defines the leading order terms to the Taylor series expansion. Insertion in Eq. 119(1 then yields 



k=0 y - 



n + k 



(z n - 1) log 



Zn - 1 



n. & 



log 



(A6) 



correct to order 0([z — 1] log|z — 1| ) . Here, 7(71, x) is the incomplete Gamma function, and can be introduced using the series 
identity 8.354.1 of Gradshteyn & Ryzhik (1980). 

Another identity required for the development of Eq. 1(28(1 pertains to the evaluation of j(n, n) for large n . For low values of 
n, say n <^ 20 , one can use 8.352.1 of Gradshteyn & Ryzhik (1980), namely 



l(n, n) = j_ _ r(n, n) = l _ e -n 



n u 

k\ 

k=0 



(A7) 



as an efficient series evaluation for either the 7(71, n) or T(n, n) incomplete Gamma functions. For larger n , one can appeal to 
the identity T(n + 1, n) = nT(n, n) + n n e~ n , together with Stirling's expansion T(n) w e~ n ri n ~ x l' 2, \/2n and identity 6.5.25 of 
Abramowitz & Stegun (1965), namely T(n + 1, n) w e~ n n n y/nir /2 , to establish the results 



7(71, n) _ 1 _ 



r(n, n) 



lim - ~ ■ ~ =: —— lim -_ 

ra^oo r(n) 2 rwoo T(n) 



(A8) 



Retaining the next order contributions to j(n, n) and r(n, n) then yields corrections of approximately ±l/(3v / 2mr) to the two 
ratios in Eq. 
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